'H-5 

INTER-ANNUAL  VARIABILITY  OF  ACOUSTIC  RAY 
TRAVEL  TIMES  IN  THE  NORTHEAST  PACIFIC 

by 

John  Alan  Furgerson 

B.S.  Elect.  Eng.,  University  of  Southern  California  (1987) 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  OCEAN  ENGINEERING 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

and  the 

WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION 

September  1990 

(c)  John  Alan  Furgerson,  1990 

The  author  hereby  grants  to  MIT,  WHOI,  and  the  U.S.  Government  permission  to 
reproduce  and  to  distribute  copies  of  this  thesis  document  in  whole  or  in  part. 


INTER-ANNUAL  VARIABILITY  OF  ACOUSTIC  RAY  TRAVEL 
TIMES  IN  THE  NORTHEAST  PACIFIC 

by 

John  Alan  Furgerson 


Submitted  to  the  Massachusetts  Institute  of  Technology/ 

Woods  Hole  Oceanographic  Institution 

Joint  Program  in  Oceanographic  Engineering 

on  August  10,  1990,  in  partial  fulfillment  of  the 

requirements  for  the  degree  of 

Master  of  Science  in  Ocean  Engineering 


Abstract 

An  acoustic  tomography  experiment  consisting  of  a  source  near  Hawaii  and  seven  receivers 
along  the  west  coast  of  North  America  was  conducted  from  November  1987  to  May  1988 
and  from  February  1989  to  July  1989.  In  this  thesis,  the  acoustic  ray  travel  times  are 
analyzed  in  order  to  investigate  inter-annual  basin-scale  thermal  variability.  These  thermal 
fluctuations  may  help  detect  any  greenhouse  warming  and  greater  understanding  of  them 
will  increase  knowledge  of  ocean-atmosphere  interactions  which  affect  weather  and  climate. 
A  discussion  of  the  program  for  finding  the  travel  times  is  included  along  with  a  comparison 
of  two  methods  of  measuring  travel  times. 
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Chapter  1 


Introduction 


1.1      Acoustic  Tomography 

Ocean  acoustic  tomography  was  first  proposed  as  a  means  of  measuring  mesoscale  (about 
100  km)  processes  by  Walter  Munk  and  Carl  Wunsch  (1979).  The  technique  involves 
measuring  the  acoustic  travel  times  between  a  source  and  receiver  and  consequently,  from 
these  measurements,  inferring  the  sound  speed  and  currents  in  between  the  instruments.  A 
1981  acoustic  tomography  experiment  (Cornuelle  et  al.,  1982;  1985)  successfully  mapped 
a  300  km  by  300  km  area  section  near  Bermuda  for  sound  speed  anomalies.  A  1981 
experiment  (Spiesberger  et  al.,  1983)  measured  Gulf  Stream  meanders. 

Acoustic  tomography  exploits  the  sound  speed  properties  of  the  ocean.  The  sound 
speed  depends  upon  the  salinity,  pressure,  and  temperature.  When  a  ray  passes  through 
an  interface  between  waters  of  different  properties  it  bend  towards  the  minimum  sound 
speed  (see  FIGURE  1-1)  according  to  Snell's  Law, 
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Figure  1-1:  Ray  traveling  through  interface  and  bending  toward  the  smallest  sound  speed 
according  to  Snell's  Law. 

where  #,■  is  the  angle  the  ray  makes  with  the  horizontal  (grazing  angle)  and  c,  is  the  sound 
speed  in  the  layer. 

In  a  typical  mid-latitude  portion  of  the  ocean,  the  speed  of  sound  decreases  with 
distance  from  the  surface  due  to  a  temperature  drop.  Below  a  certain  depth,  the  temper- 
ature is  relatively  constant  and  the  increase  in  pressure  with  depth  causes  the  sound  speed 
to  increase.  This  forms  an  axis  of  minimum  sound  speed  (sound  channel  axis)  towards 
which  the  rays  will  be  refracted.  A  sound  speed  profile  for  the  Atlantic  Ocean  is  given 
in  FIGURE  1-2  (right  side).  The  rays  will  constantly  curve  due  to  the  changing  sound 
speed.  The  paths  of  several  rays  are  drawn  in  FIGURE  1-2  (left  side)  emanating  from 
a  source  at  a  depth  of  1.2  km.  If  a  receiver  were  located  on  the  sound  channel  axis  80 
km  away,  only  certain  rays  (called  eigen-rays)  would  reach  it.  Rays  which  do  not  interact 
with  the  air-sea  interface  or  the  bottom  will  lose  very  little  energy  and  can  be  detected 
at  tremendous  distances. 
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Figure  1-2:  A  typical  sound  speed  profile  is  given  on  the  right  with  some  ray  paths  drawn 
on  the  left.  Only  certain  rays,  eigen-rays,  will  reach  the  receiver.  (From  Ewing  and 
Worzel,  1948.) 

The  advantages  acoustic  tomography  offers  over  traditional  ocean  observation  tech- 
niques include: 

1.  As  more  sources  (S)  and  receivers  (R)  are  added  to  the  tomographic  array  the 
amount  of  information  increases  by  R  X  S  instead  of  R  +  S  as  with  point  measure- 
ments. 

2.  Acoustic  tomography  comes  close  to  providing  real  time  monitoring  of  the  ocean 
since  sound  travels  at  approximately  1.5  km  per  second. 


1.2      Motivation  for  Experiment 

As  human  activities  increase  the  carbon  dioxide  (CO2)  and  other  greenhouse  gasses  in 
the  atmosphere,  most  climate  models  predict  an  increase  in  the  earth's  temperature  due 
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to  greenhouse  warming.  The  temperature  increase  would  not  be  felt  uniformly  across  the 
globe  (Stouffer  et  al.,  1989).  In  1983,  Spiesberger  began  an  acoustic  tomography  experi- 
ment for  the  primary  purpose  of  monitoring  global  warming  in  the  ocean  (Spiesberger  et 
al.,  1990a,b). 

A  problem  with  the  global  climate  models  is  the  lack  of  knowledge  of  the  ocean's 
role  in  affecting  climate.  The  upper  two  and  a  half  meters  of  the  ocean  has  the  same 
thermal  capacity  as  the  entire  atmosphere  (Gill,  1982)  and  fluxes  of  heat  transported  by 
the  ocean  and  atmosphere  are  about  equal  (Vonder  Haar  and  Oort,  1976).  To  confirm  the 
accuracy  of  the  models,  a.  reliable  method  of  measuring  the  ocean's  temperature  on  the 
same  spatial  scale  as  the  greenhouse  warming  model  predictions  is  necessary  (Barnett, 
1990).  The  model  of  Stouffer  et  al.  (1989)  suggests  a  spatial  structure  on  the  order  of 
1000  to  10,000  km. 

The  usual  method  of  determining  the  heat  content  of  the  ocean  involves  taking  point 
measurements  at  hydrographic  stations  or  by  dropping  expendable  and  air-expendable 
bathythermographs  (XBT's  and  AXBT's).  Point  measurements  of  temperature  are  con- 
taminated by  variability  on  the  meso  (~  100  km,  ~  30  days)  and  fine  (~  10  —  1000  m, 
~  0.3  —  24  hr)  scales  (Spiesberger  et  al.,  1990a, b).  Random  and  systematic  errors  are 
present  (Hanawa  and  Yoritaka,  1987)  along  with  nonuniformities  between  XBT's  from 
different  production  runs  ( Wyrtki  and  Uhrich,  1982).  Therefore,  monthly  average  temper- 
ature anomalies  cannot  be  measured  from  hydrographic  data  (Spiesberger  et  al.,  1990a,b; 
Wyrtki  and  Uhrich,  1982). 

Acoustic  tomography  is  a  powerful  tool  for  monitoring  large-scale  variability  (Spies- 
berger et  al.,  1990a, b).    Using  a  few  acoustic  instruments,  the  large-scale  thermal  fluc- 
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tuations  in  the  ocean  can  be  measured  which  may  improve  understanding  of  the  ocean's 
affect  on  weather,  climate,  ocean  circulation,  and  greenhouse  warming  (Spiesberger  et  al., 
1990a,b). 
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Chapter  2 


The  Experiment 


The  experiment  used  one  source  and  seven  receivers.  The  location  of  the  equipment  can 
be  seen  in  FIGURE  2-l(upper  left).  The  acoustic  source  was  located  on  the  Oahu  slope 
off  of  Kaneohe  Bay  in  183  m  of  water  and  the  receivers  were  bottom-mounted  at  differing 
ocean  depths. 

Source  1  (Si)  transmitted  to  all  seven  receivers  from  21  November  1987  to  9  May  1988. 
It  was  then  repaired  and  redeployed  approximately  20  m  from  its  original  position.  It 
then  functioned  from  2  February  1989  to  5  July  1989  with  a  shutdown  period  of  ten 
days  centered  around  18  April  and  one  of  five  days  around  3  May.  The  distance  from 
Si  to  the  receivers  ranged  from  2000  to  4000  km.  In  general,  data  was  obtained  from  all 
receivers  during  the  period  in  which  Si  transmitted  except  for  short  breaks  due  to  power 
interruptions  and  data  recording  failures.  No  data  is  available  in  1989  for  Rl  and  R2. 
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Figure  2-1:  Map  of  the  experiment  (upper  left)  and  the  travel  times  at  each  source-receiver 
pair  (SRP).  Travel  time  changes  (reckoned  to  the  first  1987  datum)  are  on  the  left  vertical 
axis  and  the  approximate  temperature  changes  are  given  on  the  right  vertical  axis.  The 
thin  lines  are  for  1987-88  data  and  the  thick  lines  for  1989  data. 
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Chapter  3 


Data  Processing 


3.1      Signal  Processing 

The  signal  generator  for  the  source  produces  a  pulse-like  signal.  The  signal  generator  is 
located  on  land  so  that  the  signal  must  be  amplified  before  travelling  through  a  five  mile 
length  of  cable  to  the  acoustic  source.  The  multi-path  signal  from  Si  was  recorded  every 
18  minutes  so  that  there  are  about  76  records  per  day.  The  clocks  for  the  signal  generator 
and  receivers  are  maintained  by  rubidium  oscillators  and  a  GOES  satellite  receiver. 

Since  the  source  cannot  transmit  a  pulse  powerful  enough  to  be  received  above  the 
background  noise  level,  a  phase- modulated  signal  is  sent  out  instead.  Si  used  a  133 
Hz  carrier  frequency  phase-modulated  every  8  cycles  by  a  255  digit  linear  maximal  shift 
register  sequence.  Each  digit  has  a  length  of  (8/133)  =  60  ms.  The  sequence  has  a 
period  of  255(8/133)  =  15.338  s.  By  correlating  the  received  signal  with  a  replica,  the 
energy  in  a  255  digit  period  is  compressed  into  the  60  ms  single  digit  period  increasing  the 
signal- to- noise  ratio  (SNR.)  by  101og(255)  =  24  decibels  (dB).  Four  consecutive  periods 
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are  coherently  averaged  to  increase  the  SNR  by  10 log 4  =  6  dB.  Next,  a  running  average 
of  four  complex  samples  gives  another  6  dB  gain  for  a  total  increase  in  SNR  of  about  36 
dB  (Metzger,  1983;  Bushong,  1987;  Spiesberger  et  al.,  1989). 

3.2      Stability  of  Signal 

Adjacent  individual  amplitude  records  for  yearday  153  1988  of  Source  1  to  Receiver  5 
(S1R5)  data  are  shown  in  Figure  3-1  (bottom).  The  adjacent  records  are  not  alike  due  to 
internal-wave  induced  sound  speed  fluctuations  and  interactions  with  the  bottom  (Spies- 
berger, personal  communication).  The  incoherent  averaging  of  records  over  a  given  day 
reduces  the  variability  of  waves  with  periods  much  less  than  the  given  period  (Spiesberger 
et  al.,  1980).  The  incoherent  average  is  defined  by, 

aN  =  [T^-Ell^^]||2]1/2,  (3-1) 

where  the  mth  complex  demodulate  of  the  rth  record  is  d[m,r]  and  Nrec  is  the  number  of 
records  being  averaged  (Spiesberger  et  al.,  1989).  The  incoherent  average  of  three  adjacent 
records  for  hours  one  to  six  are  given  in  FIGURE  3-l(top)  and  four  hour  incoherent 
averages  of  about  13  records  each  in  FIGURE  3-2  (bottom). 

Next  are  the  daily  (~  76  record)  incoherent  averages  for  six  consecutive  days 
(FIGURE  3-2,  top).  The  dashed  lines  indicate  the  six  rays  that  are  being  tracked  for 
this  source- receiver  pair  (SRP).  For  most  of  the  SRP's  the  daily  averages  were  used  to 
track  the  travel  times  of  the  rays.  For  S1R3,  S1R4  and  S1R7  1989  data,  two  day  averages 
(~  150  records)  were  used  to  increased  the  SNR  at  the  expense  of  halving  the  amount 
of  data.   The  final  stability  graph  (FIGURE  3-3)  for  S1R5  data  contains  daily  averages 
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Figure  3-1:  Six  adjacent  amplitude  records  (left)  and  six  incoherently  averaged  hourly 
records  (right),  consisting  of  the  three  records  which  occurred  in  the  hour,  for  S1R5 
yearday  153  1988  data.  The  records  show  little  resemblance  to  each  other  due  to  internal 
waves  perturbing  the  sound  speed  field  and  bottom  interactions.  The  background  noise 
increases  due  to  various  surface  phenomena  including  storms  (effects  last  ~  3  days)  and 
ships  (~  1  hour)  (Spiesberger  et  a!.,  1980). 


17 


^wW*^^^ 


Q 
3 


Q_ 
< 


[I'WT  HRS.    16-20 


HRS.    12-16 
HRS.   8-12 


1 1 1 ' 1 


5  6  7  8  9  10 

TRAVEL  TIME   (SEC) 


~^-^*^ 


i   i  i  1 1 
i   i 


/y^v^/w^MVi     i  iV  Vx^tf 


t 1 r — T 

5  6 


~i « r 

8  9 


TRAVEL   TIME    (SEC) 


10 


Figure  3-2:  Four  hour  (left)  and  daily  (right)  incoherent  averages  consisting  of  about 
13  and  76  records  respectively  for  S1R5  1988  data.  As  the  time  length  of  the  averages 
increase  past  the  internal  wave  periods,  the  variability  of  the  averages  decreases.  Dashed 
lines  indicate  stable  arrivals. 
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Figure  3-3:   Daily  averaged  records  from  each  of  the  three  years  of  S1R5  data.    Dashed 
lines  are  drawn  to  illustrate  the  stable  arrivals. 
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for  days  in  each  of  the  three  years.  The  daily  averages  exhibit  stability  over  multi-year 
periods. 

3.3  Ray  Travel  Time  Measurement 

A  peak  finding  program,  PCFIND,  was  used  to  measure  the  travel  times  of  the  acoustic 
rays.  Some  consider  the  peak  amplitude  to  be  the  measure  of  a  ray's  arrival  while  others 
consider  the  centroid  determined  about  the  peak  as  the  better  measure.  PCFIND  will 
find  both.  It  first  finds  the  peaks  in  the  incoherent  average  files.  A  window  is  then 
set  up  within  which  the  average  amplitude  is  determined.  PCFIND  slides  this  window 
around  the  peak  in  order  to  find  the  maximum  average  amplitude.  Within  this  window 
the  centroid  is  found.  The  travel  times  of  the  peaks  and  centroids  are  written  to  separate 
files.  Details  of  PCFIND  are  in  Appendix  A.  A  comparison  of  peak  and  centroid  data  is 
in  Appendix  B. 

3.4  Reduction  of  Data 

The  peak  travel  times  for  S1R3  are  graphed  in  Figure  3-4  including  1983  data  from  a 
previous  experiment.  To  the  first  order,  the  travel  time  changes  are  the  same  and  therefore 
the  changes  were  averaged  for  each  day.  First,  the  ray  with  the  highest  number  of  "good" 
data  points  was  chosen  as  the  reference  ray.  Then  the  average  offset  of  the  other  rays 
from  this  reference  ray  is  determined  by, 

1    N 
Toffc  =  jj  H  Tic  -  Ik,  (3.2) 
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Figure  3-4:  Peak  travel  times  of  the  five  rays  tracked  for  S1R3  including  1983  data  from 
a  previous  experiment.  The  long  breaks  in  the  data  are  due  to  source  shutdowns.  Small 
gaps  are  due  to  short  term  equipment  breakdowns  or  an  inability  to  track  the  rays.  The 
travel  times  decrease  in  1988  as  the  water  begins  to  warm  after  a  February  cold  peak.  In 
1989,  the  travel  times  increase  from  January  to  June  indicating  a  temperature  decrease. 

where  c  refers  to  the  particular  ray,  r  is  the  reference  ray,  N  is  the  number  of  common 
data  points  between  rays  c  and  r,  and  T  is  travel  time  reckoned  to  an  arbitrary  zero.  The 
average  travel  time  for  each  day  j  is  then  found  from, 

NRAYS 


T, 
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vej  "  NRA\  . 
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y^Zn     £     Tcj-Tofh,    c^r,  (3.3) 


where  NRAYS  is  the  number  of  tracked  rays.  The  standard  deviations  of  the  offsets  were 
also  obtained. 
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Chapter  4 


Results  and  Conclusions 


4.1      Travel  Times  for  Each  Source-Receiver  Pair 

The  collected  data  for  each  receiver  is  depicted  in  Figure  2-1.  In  the  upper  left  corner  of 
the  figure  is  a  map  of  the  experiment. 

For  1987-88  data,  receivers  R3,  R6,  and  R7  show  a  similar  signal.  The  waters  cool 
in  the  winter,  increasing  the  travel  times,  and  begin  to  warm  in  the  spring,  decreasing 
the  travel  times.  Rl  and  R5  in  1987-88  display  some  of  the  warming  trend  but  none  of 
the  cooling  in  the  winter.  The  travel  times  of  R2  change  very  little.  The  signal  for  R4 
is  quite  different  from  the  others.  A  slight  cooling  trend  is  seen  in  to  June  1988  when 
the  water  temperature  begins  to  increase.  In  the  beginning  of  March,  however,  a  cooling 
trend  begins  which  lasts  until  the  data  ends  in  late  April. 

All  of  the  1989  data  starts  with  a  cooling  trend  lasting  in  to  April  followed  by  a 
warming  trend.  The  1989  temperatures  become  cooler  than  the  1988  data  in  February  or 
March  for  the  SRP's  where  the  plots  cross.  R3  shows  no  clear  trend  after  mid-April.  115 
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Figure  4-1:  Spatial  average  of  all  SRP's 


A 


displays  a  cooling  trend  in  late  June. 

A  spatial  average  of  the  travel  times  of  all  SRP's  is  shown  in  FIGURE  4-1.  The  graph 
indicates  a  slight  temperature  decrease  from  November  1987  to  the  end  of  January  1988. 
The  northeast  Pacific  then  warms  from  February  to  May  1988.  As  the  1989  data  begins, 
the  waters  are  warmer  in  February  than  they  were  one  year  earlier.  A  cooling  trend  lasts 
from  February  to  April  1989  followed  by  a  slight  warming  trend.  The  overall  result  for 
the  northeast  Pacific  basin  in  1989  is  cooler  water  temperatures  in  June  than  in  February. 

4.2      Sea  Surface  Temperature  Data 

The  travel  time  data  are  compared  to  the  sea  surface  temperature  (SST)  monthly  means 
from  the  Oceanogruphic  Monthly  Summaries  (1987-  1989).  The  data  from  the  Oceano- 
graphic  Monthly  Summaries  is  a  blend  of  in  situ  data  used  to  establish  "benchmark" 
temperature  values  and  satellite  data  used  to  determine  the  shape  of  the  contour  lines 
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Figure  4-2:  Monthly  sea  surface  temperature  average  for  the  great  circle  path  between 
Source  1  and  Receiver  5. 

between  the  points.  The  in  situ  data  gives  the  bulk  temperature  of  the  top  layer  of  water 
(on  the  order  of  meters)  while  the  satellite  data  is  indicative  of  the  skin  temperature  of 
the  top  millimeter  of  water  (Reynolds,  1988).  The  satellite  data  must  be  adjusted  so  that 
it  coincides  with  the  in  situ  data. 

A  monthly  mean  SST  temperature  between  the  source  and  a  particular  receiver  was 
obtained  by  integrating  along  the  great  circle  path.  The  spatially  averaged  monthly  SST 
mean  for  S1R5  is  plotted  in  Figure  4-2  and  is  to  be  compared  with  the  travel  time  data 
in  Figure  2-1.  No  correlation  was  observed  for  S1R5  nor  for  any  of  the  other  SRP's. 

Since  the  acoustic  rays  sample  a  wider  portion  of  the  water  column  than  the  SST  data, 
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they  average  the  temperature  of  a  greater  volume  of  water.  The  tomographic  data  is  also 
less  dependent  on  seasonal  temperature  fluctuations  and  therefore  makes  a  better  measure 
of  the  inter-annual  changes  in  heat  content  of  the  oceans  (Spiesberger  et  al.,  1990a,b). 

4.3      El  Nino  Effects 

Acoustic  tomography  can  provide  extensive  temperature  maps  of  the  oceans  which  may 
provide  a  better  understanding  of  events  which  affect  weather  and  climate  such  as  Los 
Nihos.  The  term  El  Nino  originally  referred  to  the  anomalous  warming  along  the  Peruvian 
coast  which  starts  each  year  around  Christmas.  El  Nino  currently  refers  to  anomalous 
warming  centered  along  the  equator  in  the  central  and  eastern  Pacific  as  part  of  an  El 
Nino-Southern  Oscillation  (ENSO)  event.  Normally,  the  high  temperatures  in  the  central 
and  eastern  equatorial  Pacific  last  for  10-15  months  and  are  followed  by  an  anomalous  cold 
period  named  La  Nina  (Philander,  1990).  Usually,  Los  Nihos  occur  at  3-4  year  intervals 
with  stronger  events  at  least  6-7  years  apart  (Enfield,  1989).  Through  atmospheric  and 
ocean  teleconnections,  Los  Ninos  affect  areas  far  removed  from  the  equatorial  Pacific 
(Quinn  et  al.,  1987;  Philander,  1990). 

The  effects  of  El  Nino  reach  the  northeast  Pacific  basin  in  two  ways.  The  atmospheric 
link  is  provided  by  Hadley  cells  which  circulate  warm  air  rising  along  the  equator  to 
atmospheric  downdrafts  in  the  mid-latitudes.  The  ocean  link  takes  the  form  of  coastally- 
trapped  Kelvin  waves  which  move  poleward  at  about  200  cm/s  (Johnson  and  O'Brien, 
1990).  Trenberth  (1989)  credits  the  1986-1987  El  Nino  as  a  major  cause  of  the  1988 
North  American  drought.   The  cooling  of  a  La  Nina  also  affects  large-scale  atmospheric 


25 


circulation  (Bradley  et  al.,  1987). 

It  is  too  early  for  any  in-depth  analysis  of  the  northeast  Pacific  response  to  the  1986-87 
El  Nino.  Studies  of  earlier  Los  Nifios  (Emery  and  Hamilton,  1985;  Johnson  and  O'Brien, 
1990)  show  that  the  response  to  events  of  equal  intensity  can  be  quite  different.  Emery 
and  Hamilton  show  anomalous  warming  throughout  the  northeast  Pacific  nine  months 
after  the  peak  equatorial  warming  of  the  strong  1957-58  El  Nino.  But,  nine  months  after 
the  peak  equatorial  warming  of  the  strong  1972-73  El  Nino,  there  is  anomalous  cooling 
in  a  900  km  band  along  the  entire  North  American  coast  and  between  20°  N  and  30° 
N  latitude  (Emery  and  Hamilton,  1985).  After  the  moderate  intensity  event  in  1976- 
77,  there  were  above  normal  temperatures  throughout  the  tomographic  area  (Emery  and 
Hamilton,  1985). 

The  El  Nino  of  1986-87  may  have  affected  the  acoustic  travel  times.  This  event  was 
classified  as  moderate  in  intensity  by  Quinn  et  al.  (1987).  The  eastern  equatorial  Pacific 
waters  reached  a  peak  of  about  1.8°C  above  normal  in  August  1987  (Trenberth  et  al., 
1988).  Applying  a  nine  month  delay  in  the  northeast  Pacific  response  to  the  August  1987 
peak  warming,  one  would  predict  a  warming  or  cooling  peak  during  May  1988  in  the 
northeast  Pacific.  The  data  (FIGURE  4-1)  shows  warming  up  to  the  predicted  peak  in 
May  1988.  Future  work  is  required  to  determine  if  the  acoustic  travel  times  are  responding 
to  the  1986-1987  El  Nino  including  a  comparison  to  hydrographic  data. 

On  16  February  1990,  the  National  Weather  Service's  Climate  Analysis  Center  released 
an  ENSO  advisory  due  to  a  rise  in  the  SST  in  the  central  equatorial  Pacific  and  a  weakening 
of  the  equatorial  trade  winds  (Monastersky,  1990).  The  placement  of  the  receivers  may  be 
ideal  for  detecting  a  northward  moving  Kelvin  wave.  Receivers  R.6  and  R7  should  prove 
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useful  in  determining  whether  anomalous  warming  or  cooling  caused  by  the  El  Nino  is 
restricted  to  a  band  along  the  coast  or  if  it  is  basin-wide. 

4.4      Conclusions 

Acoustic  tomography  can  measure  the  thermal  properties  of  the  ocean.  The  temperature 
data  this  experiment  provides  is  uncorrelated  with  SST  data  obtained  by  standard  meth- 
ods. Since  acoustic  tomography  averages  the  ocean's  parameters  over  the  entire  source 
to  receiver  path,  it  is  a  better  indication  of  the  in  situ  heat  content  of  the  ocean  than 
satellite  data. 

The  averaging  properties  of  acoustic  tomography  should  prove  useful  in  the  measure- 
ment of  any  possible  greenhouse  warming.  The  northeast  Pacific  experiment  shows  the 
feasibility  of  detecting  inter-annual  basin-wide  temperature  changes  in  the  ocean  (Spies- 
berger  et  al.,  1990a,b).  To  detect  greenhouse  warming,  a  tomographic  system  encompass- 
ing all  of  the  world's  oceans  is  needed. 

The  1987-1989  experiment  may  indicate  part  of  the  northeast  Pacific  response  to  the 
1986-1987  El  Nino.  Acoustic  tomography  may  provide  insight  on  the  remote  effects  of 
an  El  Nino  and  in  turn  aid  in  the  understanding  of  the  inter-relationships  between  the 
atmosphere  and  oceans. 
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Appendix  A 


The  Peak  Finding  Program 


The  records  for  a  given  SRP  must  be  compared  to  each  other  visually  to  figure  out  how 
many  rays  are  stable  from  day  to  day.  The  record  with  the  most  representative  peak 
profile  is  chosen  to  be  the  replica.  The  intensity  records  are  processed  sequentially  by  the 
peak  finding  program,  PCFIND. 

The  parameters  which  PCFIND  asks  for  are  as  follows: 

NAMEREP  the  name  of  the  file  to  be  used  as  the  replica 

NRAYS  the  number  of  rays  to  search  for  in  each  record 

TMEST(NRAYS)   the  estimated  travel  time  of  each  ray  in  the  replica 

SNRLIM  ±   the  interval  over  which  an  average  noise  figure  is  determined 

MINSNR  the  minimum  signal-to-noise  ratio  of  a  peak  to  make  it  an  acceptable  data 
point 

XCORDIS  the  positive  or  negative  maximum  lag  allowed  in  the  cross-correlation  (used 


to  minimize  computer  processing  time) 

PKWIND  the  width  of  the  peak  search  window 

CENTWID  width  of  the  window  for  finding  the  maximum  average  intensity  which  is 
used  to  find  the  centroid 

CENTDIS  the  distance  to  each  side  of  the  MAX  PEAK  that  the  centroid  find  window 
is  slid  in  search  of  the  maximum  average  intensity 

The  program  first  cross-correlates  the  replica  with  a  given  record.  The  resulting  lag 
is  added  to  the  estimated  ray  travel  times  and  a  window  of  width  PKWIND,  centered  at 
each  of  these  adjusted  times,  is  set  up  to  search  for  the  intensity  peak  for  each  ray.  If  the 
maximum  intensity,  called  MAX  PEAK,  occurs  at  the  edge  of  the  window,  it  is  compared 
to  the  next  point  over  to  ensure  that  it  is  of  greater  intensity  and  indeed  a  peak.  Next, 
the  average  noise  intensity  is  determined  over  the  selected  noise  interval  and  the  peak 
intensity  is  compared  to  it.  If  the  SNR  is  less  than  the  MINSNR  or  the  MAX  PEAK  is 
not  an  actual  peak,  a  negative  number  is  recorded  in  an  output  file  for  the  travel  time 
instead  of  the  actual  time. 

The  location  of  the  MAX  PEAKS  are  the  basis  for  finding  the  centroid  travel  times. 
A  centroid  search  window  is  set  up  of  width  CENTWID.  The  center  of  this  window  slides 
a  distance  CENTDIS  to  each  side  of  the  MAX  PEAK.  At  each  data  point  that  the  center 
occupies  the  average  intensity  within  the  window  is  determined.  The  location  of  the 
centroid  search  window  center  where  the  maximum  average  intensity  was  found  is  used  as 
the  balance  point  in  calculating  the  centroid  within  the  window.  These  are  called  MAX 
CENTROIDS  and  are  recorded  in  a  separate  file  from  the  MAX  PEAKS. 
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Another  program  filters  out  travel  times  for  each  ray  which  are  unreasonably  different 
in  value  from  the  other  data  for  that  ray.  The  program  compares  each  ray's  travel  time 
to  a  running  average  of  the  previous  five  days  of  data  for  that  ray.  If  the  difference  is 
greater  than  a  chosen  value  (usually  80  msec),  the  data  point  is  discarded. 
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Appendix  B 


Comparison  of  Peak  and 


Centroid  Data 


The  intensity  record  is  corrupted  by  noise  so  that  the  exact  arrival  time  of  an  acoustic  ray 
is  unknown.  One  method  of  tracking  the  arrivals  is  to  assume  that  the  time  associated 
with  the  highest  intensity  peak  is  the  travel  time.  Since  the  peaks  are  not  fine  lines,  some 
consider  the  centroid  about  the  maximum  intensity  to  be  the  best  estimate  of  the  arrival 
time. 

To  show  the  difference  between  what  PCFIND  determines  to  be  the  MAX  PEAK  and 
the  MAX  CENTROID,  the  intensity  record  for  day  213,  hour  18  of  S1R6  data  is  given 
in  FIGURE  B-l.  In  this  figure,  rays  2  and  8  are  marked  and  shown  in  more  detail  in 
FIGURE  B-2.  The  difference  in  peak  and  centroid  travel  times  for  ray  2  is  42  msec  and  for 
ray  8,  3  msec.  The  situation  for  ray  8  where  the  MAX  CENTROID  and  PEAK  essentially 
coincide  is  the  usual  one.    For  ray  2,  there  is  more  energy  situated  near  11.32  sec  than 
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Figure  B-l:  Intensity  Record  for  S1R6,  day  213  of  1987,  hour  18.  Rays  2  and  8  will  be 
examined  to  show  how  the  peak  finding  program  determines  the  peak  and  centroid  travel 
times. 
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Figure  B-2:  In  the  top  graph,  the  travel  times  are  determined  for  ray  2.  The  program 
looks  in  the  peak  find  window  for  the  maximum  intensity.  This  is  found  at  about  11.360 
seconds  and  is  labelled  MAX  PEAK.  Then  a  centroid  find  window  is  formed  which  slides 
along  the  time  axis  a  distance  called  the  centroid  slide  width  on  each  side  of  MAX  PEAK. 
Within  this  centroid  find  window,  the  average  intensity  is  determined.  The  time  location 
of  the  centroid  find  window  middle  with  the  maximum  average  intensity  will  be  used  as 
the  balancing  point  in  the  determination  of  the  centroid.  The  centroid  of  this  maximum 
average  intensity  window  is  found  at  about  11.318  seconds  and  is  labelled  MAX  CEN- 
TROID. In  the  bottom  graph  the  travel  times  for  ray  8  are  shown.  This  is  the  more  usual 
case  where  the  MAX  PEAK  and  MAX  CENTROID  are  close  together. 
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Type  of  Data 

Number 
of  Data 

Average  Travel 
Time  (s) 

Standard 
Deviation  (ms) 

Mean  Standard 
Deviation  (ms) 

peak 

7 

11.312 

:v> 

12 

centroid 

7 

11.308 

26 

10 

Table  B.l:  Means  and  Standard  Deviations  of  S1R6  Ray  2  1988  Yearday  213 


Type  of  Data 

Number 
of  Data 

Average  Travel 
Time  (s) 

Standard 
Deviation  (ms) 

Mean  Standard 
Deviation  (ms) 

peak 

6 

11.304 

27 

11 

centroid 

6 

11.306 

27 

11 

Table  B.2:  Means  and  Standard  Deviations  with  Hour  18  Data  Removed 

11.36  sec  so  the  MAX  CENTROID  is  determined  in  the  former  area. 

The  centroid  search  window  may  have  been  allowed  to  slide  too  far  from  the  MAX 
PEAK  in  finding  the  maximum  average  intensity.  To  find  out,  the  means  and  standard 
deviations  of  ray  2  for  yearday  213  must  be  examined  (TABLE  B.l).  The  third  column 
represents  the  number  of  "good"  travel  time  data  points  out  of  the  twelve  records  for  that 
day.  The  mean  standard  deviation  is  given  by: 


MSD  =  flNDATA) 


(B.l) 


where: 
MSD  the  mean  standard  deviation 


SD 


the  standard  deviation 


NDATA     the  number  of  "good"  data  points 

TABLE  B.2  then  gives  the  mean  and  standard  deviation  for  ray  2  with  the  hour  18 
data  removed.  The  MAX  CENTROID  time  for  yearday  213  hour  18  seems  much  more 
plausible  than  the  MAX  PEAK  time  when  compared  to  the  average  travel  times  for  ray 
2  for  that  day.    Additionally,  the  mean  standard  deviation  decreases  for  the  peak  data 
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when  hour  18  data  is  removed  while  the  centroid  mean  standard  deviation  increases  which 
shows  that  the  centroid  data  is  the  better  estimate  and  PCFIND  did  not  conduct  too  wide 
of  a  search  for  the  maximum  average  intensity. 

Five  days  of  S1R6  data  were  analyzed  to  decide  if  the  MAX  CENTROID  or  MAX 
PEAK  data  is  the  better  estimator  of  the  acoustic  ray  arrival  times.  The  averages  and 
standard  deviations  for  the  S1R6  data  are  shown  in  next  three  tables.  TABLES  B.3  and  B.4 
contain  the  MAX  PEAK  and  MAX  CENTROID  data  respectively  and  TABLE  B.5  con- 
tains the  same  data  for  the  time  differences  (MAX  PEAK  -  MAX  CENTROID).  The 
standard  deviations  of  the  difference  data  are  lower  than  the  PEAK  and  CENTROID  val- 
ues which  indicates  that  the  two  types  of  data  attempt  to  measure  the  same  parameter. 
The  larger  standard  deviation  values  in  TABLE  B.5  indicate  that  for  one  of  the  particular 
ray  travel  time  measurements  for  that  day  there  was  the  situation  where  the  MAX  PEAK 
and  CENTROID  were  far  apart  as  in  the  previous  hour  18  ray  2  data.  The  average  of  the 
mean  standard  deviation  for  the  five  days  is  9.35  msec  for  the  peak  data  and  8.95  msec 
for  the  centroid  data.  The  centroid  data  would  therefore  be  the  better  measure  of  travel 
time  but  only  because  of  the  5.5%  difference  in  mean  standard  deviations.  Since  only  the 
trends  in  the  travel  time  changes  and  not  the  absolute  travel  times  are  examined,  MAX 
PEAK  data  has  been  used  for  all  the  graphs  since  it  is  easier  to  visualize  the  intensity 
peak  as  the  mark  of  an  acoustic  ray's  arrival. 
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YRDAY 

RAY  # 

NUMBER  OF 
DATA 

213 

1 

5 

213 

2 

7 

213 

3 

.  7 

213 

4 

a 

213 

5 

7 

213 

6 

a 

213 

7 

7 

213 

0 

5 

217 

1 

5 

217 

2 

6 

217 

3 

6 

217 

4 

5 

217 

5 

a 

217 

6 

7 

217 

7 

6 

217 

0 

6 

221 
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9 

221 

2 

10 

221 

3 

9 

221 

4 

11 

221 

s 

0 

221 

6 

V 

221 

7 

a 

221 

B 

10 

225 

1 

6 

225 

2 

6 

225 

3 

e 

225 

4 

6 

225 

5 

7 

225 

6 

6 

225 

7 

S 

225 

G 

a 

229 

1 

3 

229 

2 

4 

229 

3 

6 

229 

4 

V 

229 

5 

<? 

229 

6 

12 

229 

7 

9 

229 

0 

a 

AVG.  TIME  (S) 

.1105660E+02 
.1131157E+02 
.120180CE+02 
.1227337E+02 
.1288514E+02 
.1310812E+02 
.1362486E+02 
.1383760E+02 
..  1106160E+02 
.1131633E+02 
-1201683E+02 
. 1227740E+02 
.1288337E+02 
.1307471E+02 
-1360650E+02 
.1382467E+02 
.1103733E+02 
.1127800E+02 
•1198700E+02 
. 1224364E+02 
.1283950E+02 
. 1309033E+02 
.1357887E+02 
.1380560E+02 
.1100733E+02 
.1126550E+02 
.1195775E+02 
.1222483E+02 
.1284000E+02 
. 1305450E+02 
•1358150E+02 
. 1380762E+02 
•1101667E+02 
.1125375E+02 
-1195567E+02 
. 1219056E+02 
.1282622E+02 
.1306258E+02 
. 1356400E+02 
-1377637E+02 


STANDARD 
DEVIATION  (S) 
. 1774933E-01 
.3202157E-01 
.1801585E-01 
.1163990E-01 
.3937890E-01 
.3785970E-01 
. 23883 16E-01 
.2257069E-01 
.1657228E-01 
.1769807E-01 
•1766732E-01 
.2668786E-01 
.2949956E-01 
.3012197E-01 
•2176961E-01 
.1407929E-01 
.1519518E-01 
.2900353E-01 
.3047766E-01 
•2010084E-01 
.2961860E-01 
.2456738E-01 
.3616437E-01 
.2896625E-01 
.1457937E-01 
.1434976E-01 
.2774777E-01 
. 26302 15E-01 
.2096930E-01 
.2690580E-01 
.2147081E-01 
•1299946E-01 
•2003886E-01 
.4865908E-01 
.2242501E-01 
.4202406E-01 
.1045390E-01 
.3208445E-01 
.1335841E-01 
.4185682E-01 


MEAN  STANDARD 
DEVIATION  (S) 
.  7937743E-02 
.1210302E-01 
.6809351E-02 
•4115327E-02 
.1488382E-01 
. 1338543E-01 
•9026988E-02 
.  1009392E-01 
.7411350E-02 
.7225205E-02 
.7212652E-02 
.1193517E-01 
•1042967E-01 
•1138504E-01 
.8887408E-02 
.5747848E-02 
.  5065060E-02 
.9171725E-02 
•1015922E-01 
.  6060630E-02 
.1047176E-01 
.8189127E-02 
.1278604E-01 
.9159934E-02 
.5952004E-02 
.5858266E-02 
.9810320E-02 
.1073781E-01 
.7925649E-02 
.1098425E-01 
.7591077E-02 
.4596004E-02 
.1156944E-01 
.2432954E-01 
.9154974E-02 
.  1400802E-01 
.3484634E-02 
.9261985E-02 
.4452802E-02 
.1479862E-01 


Table  B.3:  Peak  Means  and  Standard  Deviations 
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4 
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5 
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6 
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7 

9 
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AV6.  TIME  (S) 

-1105800E+02 
.1130786E+02 
.1201686E+02 
.1227325E+02 
. 1288943E+02 
.1310787E+02 
.1362186E+02 
.1383800E+02 
.1106600E+02 
.1131867E+02 
-1202317E+02 
. 1227860E+02 
. 1288362E+02 
.1307971E+02 
. 1360250E+02 
. 1382483E+02 
.1103811E+02 
.1128160E+02 
.1199156E+02 
. 1224609E+02 
.1283638E+02 
•1309033E+02 
.1357350E+02 
. 1380470E+02 
•1100500E+02 
.1126667E+02 
.1195150E+02 
. 122226 7E+02 
. 1283929E+02 
. 1305717E+02 
•1358700E+02 
.1381025E+02 
.1100233E+02 
.1125100E+02 
•1195717E+02 
.1219722E+02 
.1282167E+02 
. 1306358E+02 
.1356478E+02 
. 1377675E+02 


STANDARD 
DEVIATION  (S) 
. 1802216E-01 
.2572620E-01 
.1730776E-01 
.2449362E-01 
.3871459E-01 
.4390460E-01 
.2440396E-01 
•2394151E-01 
•1019820E-01 
.1937057E-01 
.1975206E-01 
.2452434E-01 
•2853480E-01 
.2715423E-01 
.1161540E-01 
.1201946E-01 
.1403538E-01 
.3073500E-01 
•3005962E-01 
.1788062E-01 
.3229141E-01 
.2386558E-01 
.3477056E-01 
.2675831E-01 
.1329172E-01 
. 1387644E-01 
.2876630E-01 
.3127135E-01 
•1973917E-01 
.1628318E-01 
.2028546E-01 
.1335817E-01 
.4027969E-02 
.4903077E-01 
.2308982E-01 
•3415251E-01 
.1804918E-01 
.3148407E-01 
•1559771E-01 
.4357379E-01 


MEAN  STANDARD 
DEVIATION  (S) 
.8059755E-02 
.9723590E-02 
.6541718E-02 
.8659803E-02 
.1463274E-01 
.1552262E-01 
•9223829E-G2 
.1070697E-01 
.4560774E-02 
.7908001E-02 
.8063747E-02 
. 1096762E-01 
.1008858E-01 
.1026333E-01 
.4741968E-02 
.4906926E-02 
.4678459E-02 
.9719262E-02 
.1001987E-01 
.5391211E-02 
.1141674E-01 
.7955195E-02 
•1229325E-01 
.8461723E-02 
. 5426320E-02 
. 5665034E-02 
.1017042E-01 
.  1276648E-01 
.7460704E-02 
.6647579E-02 
.7171993E-02 
.4722827E-02 
.2325549E-02 
.2451539E-01 
•9426382E-02 
.1 13841 7E-01 
•6016394E-02 
.9088671E-02 
.5199237E-02 
.1540566E-01 


Tabic  B.4:  Centroid  Means  and  Standard  Deviations 
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4 

9 
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5 
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229 

6 

12 

229 

7 
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229 
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AVG.  TIME  (S) 

-.14O0000E-O2 

-3714286E-02 

.1142857E-02 

.1249997E-03 

-.4285715E-02 

•2500001E-03 

. 3000000E-02 

-.3999999E-03 

-.4400O00E-02 

-.2333333E-02 

-.6333333E-02 

- . 1200000E-02 

-.2500000E-03 

-.5000000E-02 

.4000000E-02 

-.16&6666E-03 

-.7777777E-03 

- . 3600000E-02 

-.4555555E-02 

-.2454545E-02 

.3125000E-02 

.3880511E-10 

-5375000E-02 

.9000000E-03 

.2333333E-02 

-.1166667E-02 

.6250000E-02 

.2166666E-02 

•7142857E-03 

-.2666667E-02 

-.550OO00E-O2 

-.2625000E-02 

.1433333E-01 

.2750000E-02 

■-1500000E-02 

•.6666667E-02 

.4555555E-02 

-.1000000E-02 

•.7777778E-03 

•.3750000E-03 


STANDARD 
DEVIATION  (S) 
.6343502E-02 
•1597958E-01 
.4940464E-02 
•1650331E-01 
.2490800E-02 
.1531952E-01 
.1511858E-02 
.4586938E-02 
. 1059434E-01 
.2624670E-02 
. 1278889E-01 
.4069398E-02 
. 3699663E-02 
.1109698E-01 
.1690168E-01 
.4669643E-02 
.3705185E-02 
.6311893E-02 
•8744668E-02 
•4  599677E-02 
•1044555E-01 
.3915782E-02 
.8092245E-02 
-9417538E-02 
. 3299833E-02 
.3890872E-02 
.1503953E-01 
.1473564E-01 
•3057276E-02 
.2333571E-01 
•1013657E-01 
•3870966E-02 
•1602776E-01 
.3561953E-02 
•3403430E-02 
.1235584E-01 
.1340905E-01 
.2581989E-02 
.3614102E-02 
.4498264E-02 


MEAN  STANDARD 

DEVIATION  (S) 

.2e36901E-02 

•6039714E-02 

.1867320E-02 

.5834803E-O2 

.9414338E-03 

.5416269E-02 

.5714286E-03 

.2051341E-02 

.4737934E-02 

.1071517E-02 

.5221042E-02 

•1819890E-02 

.1308029E-02 

.4194264E-02 

.  6900082E-02 

.  1906374E-02 

.  1235062E-02 

.  1995996E-02 

.2914890E-02 

.1386855E-02 

.3693058E-02 

.  1305261E-02 

.2861041E-02 

.2978087E-02 

.1347151E-02 

. 1588442E-02 

•5317279E-02 

.6015799E-02 

•1155542E-02 

.9526766E-02 

-3583819E-02 

•1368593E-02 

•9253630E-02 

.1780976E-02 

-1389445E-02 

.4118613E-02 

.4469683E-02 

•7453561E-03 

.1204701E-02 

.1590377E-02 


Table  B.5:  Peak  -  Centroid  Means  and  Standard  Deviations 
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